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Q\ , Abstract 

The thermodynamic formalism of Ruelle, Sinai, and Bowen, in 
which chaotic properties of dynamical systems are expressed in terms 
of a free energy-type function ip(/3), is applied to a Lorentz Lattice 
Gas, as typical for diffusive systems with static disorder. In the limit 
of large system sizes, the mechanism and effects of localization on 
large clusters of scatterers in the calculation of ip(P) are elucidated 
and supported by strong numerical evidence. Moreover it clarifies 
and illustrates a previous theoretical analysis (l| of this localization 
phenomenon. 
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1 Introduction 

Many nonlinear physical problems involve a complicated discrete distribu- 
tion function {pi}. As it may vary in a very irregular way with i, it would 
be attractive to replace the p^s by a smooth function containing the same 
information about the structure of the distribution. One way to do it is 
to associate with the distribution a whole set of so-called "escort distribu- 
tions" H defined as Pi = (pi) 13 /Z(/3), where Z(f3) is a partition function 
Z(0) = Sf(Pi) • The parameter (3 allows one to scan the structure of the 
initial distribution. Large (3 values enhance the most probable trajectories 
whereas negative /3's focus on the least probable trajectories (we impose that 
{PiY = if pj = so that our definitions still hold for negative /3's). By anal- 
ogy with thermodynamics where (3 would be an inverse temperature, a free 
energy — like function ip(/3) is introduced, which is related to the logarithm 
of the dynamic partition function Z(/3). 

This formalism introduced by Ruelle, Sinai and Bowen || is called ther- 
modynamic formalism. It has been successfully applied for instance to mul- 
tifractals H. 
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The present paper deals with another frame of application of this formal- 
ism, i.e. the chaotic properties of dynamical systems. For a given map (we 
will assume hyperbolicity in order to ensure good ergodicity properties), the 
phase space is partitioned into cells. Each sequence of cells explored by a 
trajectory (r Q , • • • , r t ) in t time steps is one point fl t of the dynamical phase 
space. We shall refer to it as a trajectory over t time stepsQ. 

The thermodynamic formalism will be applied to the distribution P(fl t ), 
where fit replaces the previous subscript i. The partition function reads 
now Z(fl,t) = J2n t P^iP-t)- Note that each fl t is counted with equal weight, 
i.e. the calculation does not require the a priori knowledge of the invariant 
measure. The free energy is nowf] defined as ip(/3) = lim^oo | lnZ(/3). In 
this context, it is called Ruelle pressure or topological pressure. It contains 
information about the dynamics of the system, which can be extracted again 
by varying the control parameter {3. For example, for an open system, the 
escape rate 7 is such that ^(l) = —7, and ip((3) vanishes when (3 is equal to 
the fractal dimension of the repeller (set of trajectories which never escape). 

It has already been found in a more general frame that some non analyt- 
icity of the Ruelle pressure may occur at certain (3 values |2|, ||, f| . They are 
called phase transitions in analogy with thermodynamics. 

In this paper, we study the influence of disorder on the structure of the 
Ruelle pressure. We have applied the thermodynamic formalism to systems 
with static disorder and have found very peculiar features ||. In the limit 
of infinite systems, and for almost all (3 values, the Ruelle pressure becomes 
completely determined by trajectories localized on rare fluctuations of the 
disorder. The global structure of the disorder — in particular the density of 
impurities — becomes irrelevant. Depending on (3, the relevant trajectories 
are not always localized onto the same type of fluctuations. The crossovers 
between the different localization regimes are characterized by non-analytical 
points in the Ruelle pressure. 

For finite but large systems, the Ruelle pressure is again analytic. How- 
ever for most /3's it is still dominated by localized trajectories. Delocalized 
trajectories become relevant only in some limited (3 regions which shrink to 
points as the system size tends to infinity. 

To be more precise, we have developed these ideas by studying a par- 
ticular model, the Lorentz Lattice Gas (LLG). A Lorentz gas consists of a 
moving light particle scattered among fixed heavy particles. It may be con- 
sidered as a simple model for diffusion, electric conductivity or flow inside a 
porous medium. In the past decades, it has been successfully used to study 
the connection between the irreversible behavior of fluids and their chaotic 
properties at a microscopic level. For example, P. Gaspard and G. Nicolis 
have shown that for an open system, the diffusion coefficient can be ex- 
pressed in terms of the Kolmogorov-Sinai entropy and the sum of the positive 
Lyapunov exponents for non escaping trajectories. The name Lorentz lattice 
gas is used when the light particle is, for reasons of simplicity, constrained 
to move on a (cubic) lattice, with scatterers located at the nodes -or sites. 

The localization processes occurring when the thermodynamic formalism 

2 Mathematically,this could also be formulated in terms of cylinder sets, i.e. the set of 
initial conditions that follow the trajectory Q t during the first t time steps. 
3 For systems with escape, Z(f3, t) vanishes exponentially. 
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is applied to such a system, have already been discussed in an extensive the- 
oretical analysis [|IJ. The present paper contains the numerical counterpart 
and demonstrates the validity of our analytical results. It clarifies the mech- 
anism of localization in large systems by explicitly calculating the topological 
pressure for specific non-random configurations in section |3|. In section |] an 
exact expression for the Hausdorff dimension of the repeller is obtained for 
an open LLC Moreover, numerical studies allow us to find a good estimate 
for the Ruelle pressure for large but finite systems. As an intermediate result, 
the distribution for the largest cluster of scatterers in random configurations 
(section ^ is estimated theoretically, and measured through numerical sim- 
ulations. 

2 Thermodynamic formalism 

First we give a more precise definition of our model. A number of iV fixed 
scatterers are randomly placed with a probability p on the L sites of a finite 
one-dimensional lattice T>, having either periodic or absorbing boundaries. 
The presence of a scatterer at site r = {1,2, ■ • • , L} is indicated by the 
boolean variable p(r), equal to 1 (0) if the site r is occupied (empty), with 
(p)=P = N/L. ' 

A light particle, moving on the lattice, is specified at time t (t = 0, 1, 2, • • •) 
by a state x t = {r t ,c it }, where its position r is a site on the lattice and 
its velocity, Cj = ±1, connects site r t to one of its nearest neighbors. Let 
(f>(x,t) = <pi(r,t) be the probability to find the moving particle in state x = 
{r,Ci}. 

The time evolution of the particle from time t to t+1 consists of a possible 
collision followed by propagation. The collision is defined by the following 
rules : if there is no scatterer at r t , then the particle moves ballistically, i.e. 
r t+ x = r t + Cj. This is expressed in terms of 0» as : 

<t>i(T + d,t + l) = <Mr,t). (1) 

If there is a scatterer at r t , i.e. if p(r 4 ) = 1, then the velocity of the particle 
is reversed with probability q (c' it = —c it ) or left unchanged with probability 
p (c' it = c it ), with p + q = 1. In the propagation step the particle moves 
over one lattice distance in the direction of its post-collision velocity, i.e. 
r f+ i = r t + c' it . Then, it hops from one site to the next one in the direction 
of its velocity Cj. Again, the corresponding evolution of <pi is given by : 

0»(r + a,t + l) = (uf>i(r, t) + 60_i(r, t) (2) 

with a = p and b = q = 1 — p (the reason for this notation will become clear 
later). More generally, we define site dependent transition probabilities 

a(r) = ap(r) + 1 — p(r) 

S(r) = bp(r). (3) 
They depend on the precise configuration of scatterers under consideration. 
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The probability <p(x, t) evolves with time according to a Chapman-Kolmogorov 
(CK) equation with site dependent transition probabilities, obtained by com- 
bining equations (P and (Q) : 

0,(r + a,t + l) = a(r)&(r, t) + S(r)0_<(r, t), (4) 

and more formally 

<j>{x,t + l) = Y,w(x\y)<l>{y,t). (5) 
y 

The transition matrix w(x\y) represents the probability to go from state 
y = {r', Cj} to state x = {r, Cj}, and is given by 

w(x\y) = 5(r- c h r') <%a(r') + 5 i(i+1 6(r') . (6) 

In the case of absorbing boundaries (open system), boundary states y = 
{r',Cj} = {1,+} and {L, — } referring to a particle entering the domain T> 
are excluded from the y-summation. This is equivalent to imposing the 
absorbing boundary conditions (ABC) 

0+(l,t) = 0_(L,t)=O. (7) 

In the case of periodic boundaries (closed system), we impose the periodic 
boundary conditions (PBC) 

&(r + L,t) = &(r,t). (8) 

The transition matrix satisfies the normalization relations 

ix*){<; { tzl < 9 > 

The inequality sign in Eq.(|9|) for open systems refers to the case where 
y = {r, Cj} denotes a state at a boundary site r with non-entering veloc- 
ity (boundary states with entering velocity do not occur). Indeed, the sum 
over x excludes states where the particle has escaped from the domain T>. 
Hence the probability for remaining inside the domain decreases when the 
particle finds itself on a boundary site. 

The special case p = 1 is called a Persistent Random Walk (PRW) where 
a(r) = p and b(r) = q for all r. Then the moving particle simply undergoes 
a random walk with correlated jumps. 

For later analysis it is convenient to also have a different representation of 
the CK-equation. It can be obtained by eliminating the probabilities 4>i(r, t) 
at non-scattering sites with the help of Eq. P). Let (k = 1, 2, • • ■ , N) be 
the position of the k-th scatterer and = r^ +1 — the free interval between 
scatterers. For ABC we define in addition R a = r\ — 1 and = L—r^. Then 
it is straightforward to show that the scattering amplitudes Ui(k, t) = 4>i(r k , t) 
(probability at scattering site r^) satisfy the closed set of equations 

U+(k + l,t + R k ) = aU + (k,t) + bU„{k,t) 
U-(k - l,t + Rk-!) = bU+{k,t) + aU-(k,t) (10) 
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with the boundary conditions 



Ui{N + k,t) = Ui{k,t) (PBC) 
U+(l,t) = U-(N,t) = Q (ABC). (11) 

Note that Eg . (|T0f) only depends on the set of random intervals {Rk \ k = 

0. 1, 2, • • • , N} (the coefficients a = p and b = q are sure variables), whereas 
the CK-equation (f|) depends on the set of random transition rates (a(r), 6(r) | 

1, 2, • • • , L}. It will appear below that the basic chaos properties can be ex- 
pressed in terms of the largest eigenvalue A and the corresponding left and 
right eigenvectors v (x) = fj(r) and u(x) = Ui(r) of the non-symmetric matrix 
w appearing in the CK-equation (|5|) : 

wu = Au 

vw = Av. (12) 

Then, any solution <f)(x,t) of the CK-equation approaches A l u(x) for t — > 

00. For our purpose it is again more convenient to deal with Eq.flHH) and 
determine only the components Ui(k) of eigenvectors at the scattering sites, 

1. e. Ui(k,t) = A t Ui{k) and the eigenvalue equation follows from (|10|), i.e. 



A Rk U+(k + 1) = aU+(k) + bU-(k) 
A^U^k - 1) = bU + (k) + aU„(k), (13) 

where the components Ui(k) satisfy the boundary conditions (|TTJ) . 

To describe the thermodynamic formalism we introduce a dynamical 
phase space consisting of all possible sequences Q t = {^1,^2, • • • ,x t }, which 
represent an allowed (i.e. in the ABC case, non-escaping from domain T>) 
of the moving particle visiting the state x T = {r T , c r } at the r-th time step. 
The probability P(p, t \x ) on Q t , given that the moving particle is in x Q at 
r = 0, is given by the multi-time distribution function 

P(n t \x ) = U t T - = 1 w(x T+1 \x T ) (14) 

on account of the CK-equation (||). 

The dynamic partition function is then introduced as a sum over state Qt 
in this dynamical phase space : 

Z(j3,t\x a ) = Y,i p (^t\x )f 



J2 K~^w{(3;x T+1 \x T ) (15) 
y 

In analogy with the methods of equilibrium statistical mechanics, there is an 
inverse temperature-like variable (3, which allows one to scan the structure 
of the probability distribution for Q. On the second line of (15) we have 
introduced the pseudo-transfer matrix, as 

w((3;x\y) = (w(x\y)f (16) 
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and w t (/3) denotes the t-th power of matrix w{(3). The largest eigenvalue 
A(/3) and corresponding left and right eigenvectors v((3,x) and u(/3,x) of 
w(/3;x\y) are defined analogously to (0), where a and b in equations @- 
(]T3|) take the values : 

a = pP; b = q f3 . (17) 

As the system is ergodic ||, Z(j3,t\x ) does not depend on the initial 
condition x Q in the long time limit, for almost all configurations of scatterers. 
We already mentioned that Z(j3, t\x Q ) vanishes for open systems as time tends 
to infinity More precisely, for large times, the sum ([15]) is dominated by the 
largest eigenvalue A(/3) of the pseudo transfer matrix w(/3; x\y), which is 
non-degenerate for ergodic systems. More explicitly, we use the spectral 
decomposition of u>*(/3), i.e. 

w\p- x\y) = v({3, x) [A{(3)f u(/3, y) + £ MP, x) [A n (/3)f u n ((3, y), (18) 

where A n < A for all n ^ 0. Thus the second term decays exponentially 
faster than the first one, and we obtain for large t : 

Z{/3,t\x ) ~ [A(/?)f ^2v((3, y)u{(3, x ). (19) 
y 

It should be noted that the partition function depends on the configuration 
of scatterers under consideration. 

The Ruelle or topological pressure ip((3,p) is defined as the infinite time 
limit of the logarithm of Z per unit time step, in a way similar to the definition 
of the free energy per particle in the canonical ensemble in the thermody- 
namic limit, i.e. 

^{(3 )P )=\im-(\nZ{f3 1 t\x )) (20) 

where (. . .) p indicates not only an ensemble average over all initial condi- 
tions, but also over all realizations of the disorder, i.e. all configurations of 
scatterers. Again, it can be expressed in terms of the largest eigenvalue of 
the matrix w{(3) as 

V-OM = <toA(/3))„ (21) 

where we have taken the infinite time limit inside the configurational average. 

For some specific (3 values, the Ruelle pressure has an explicit physical 
meaning : the positive Lyapunov exponent is A = — 0'(1); the escape rate 
for open systems is 7 = — -0(1); the Kolmogorov-Sinai entropy hxs follows 
from the generalization of Pesin's theorem and yields h^s — 0(1) — "0'(1); 
the topological entropy h T satisfies h T = 0(0); the Hausdorff dimension dn 
of the repeller (the set of trajectories that never escape) for an open system 
is the zero-point of the Ruelle pressure, i.e. ip{dn) = 0. A prime in the 
above formulas denotes a /3-derivative. 

For an open system, the transition matrix w is not stochastic, i.e. its 
largest eigenvalue is strictly smaller than one. Due to the loss of trajectories 
at each time step, the eigenvectors for (3 = 1 decay according to 

wu = exp(— j)u 

vw = exp(— j)v. (22) 



6 



In order to obtain the invariant vector 7T, the eigenvectors have to be rescaled 
at each time step, as explained in Ref. ||. A new transition matrix II is 
defined as 

U(x\y)=exp( 7 )w(x\y)'^. (23) 

This matrix is stochastic, i.e. its largest eigenvalue is equal to 1 associated 
to a left eigenvector x( x ) = 1, as can be seen from : 

£ X (x)n(x|y)=£n(x|y) = l. (24) 

X X 

The corresponding right eigenvector is the invariant vector 

-« = W < 25 > 

where it can be verified that 

J2^(x\y)Hy) = rr(x). (26) 
y 

With definition (^5|), 7r is normalized. The invariant vector 7r gives the proba- 
bility of finding the particle on a given site and with a given velocity, provided 
it has not escaped after infinite time. 



3 Chaos properties of special configurations 
3.1 Mean— field configurations 

In the subsequent subsections we develop a theoretical understanding of some 
typical properties of configurations of scatterers, which are relevant for de- 
scribing the dynamic partition function and Ruelle pressure in large systems. 

This will be done by a theoretical analysis of the largest eigenvalue of a 
number of relevant configurations. First we discuss in subsection p.l| mean 
field-type configurations, relevant for (3 close to unity, as discussed in ref- 
ill] . Here the fluctuations in the lengths of the interval between 



erence 



scatterers are small. In the next subsections we study configurations with an 
increasing number of solid clusters of scatterers (regions of density p = 1), 
separated by regions free of scatterers (voids) , as illustrated in figure [3]. This 
is helpful in order to understand the mechanism of localization. Large voids 
have the tendency to divide the system in independent subsystems with a 
higher density of scatterers. We will show that only the subsystem containing 
the largest cluster is relevant for determining the dynamic partition function. 

We start by considering periodic arrays of scatterers with a constant free 
interval Rk = R for k = 1,2, ■ ■ ■ , N which corresponds to a PRW on an 
ordered lattice with lattice distance R. The eigenvalue equation flT3"|) can be 
solved by making the ansatz Ui(k) = Aiexp(tqk). By setting the resulting 
secular determinant in (|T3|) equal to zero, one finds for the largest eigenvalue 
A : 

[A(/5)] fl = acosg+ b 2 — a 2 sin 2 q 

~ (a + 6){l-|-g 2 + ---j (q small). (27) 
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Figure 1: Three examples of configurations with N scatterers on a lattice 
of size L: (a) mean-field configuration; (b) one solid cluster; (c) two solid 
clusters of size M and M = N — M. Above (below) each line of sites, we 
indicate the labeling of sites (scatterers). 



The wave number q has to be determined from the boundary conditions 
(PBC or ABC). In the PBC case the allowed g-vectors follow from Ui(k + 
N) = Ui(k), so that q = 2nn/N with n = 0, ±1, ±2, ■ ■ -. Consequently, the 
wavenumber q = yields the largest eigenvalue, [A(/5)] i? = a + b. 

Next consider the open system, with ABC ( U + (l) = U-(N) = ). The 
eigenvector Ui(k) is a special linear combination of exp(iqk) and exp(— iqk), 
i.e. 

U + (k) = Asinq(k-l) 

U-{k) = B sin q(N-k). (28) 



To determine the allowed g-value we substitute equation (]2lf ) into the equa- 



tions (13), taking k — 1 for the first one and k = N for the second one. This 
yields : 

Advising = bB sin q(N-l) 

A R Bsinq = bA sin q(N - 1). (29) 

The ratio of the two equations yields A = ±B. As the components (|28|) of 
the largest eigenvector of the positive matrix w{f3) need to be positive (they 
represent probability densities), it follows that A = B = 1. Substituting 
expressions ( ]2"8| ) in the first equation of ( ]TB| ) with k = 2, and using equation 
(|29| ) to eliminate A R from the left hand side, we obtain a transcendental 



equation for the wavenumber q 

asinq = bsin(qN) (30) 

valid for = 1, 2, • • • with ABC. For large N, the smallest root of this 
equation is close to q ~ tt/N; so we set qN = n — 5q such that sin(gA^) = 
sin(g5) and find for small q that 5 = a/b. This yields for the ABC case the 
smallest allowed wavenumber : 

71 71 

* = W+S = N + a/b (31) 
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and the largest eigenvalue : 

[A(P)] R = (a + b){l 




71 



N + a/b / 



+ 0(iV- 4 )} 



(32) 



From the summary below equation (pi]) , all chaos quantities for periodic 
arrays of scatterers can be calculated from (|3"2"|). The mean field theory for 
the LLG, discussed in fLOf , follows from this result by setting R equal to 
the average free interval length (mean free path) R = L/N = 1/p and the 
resulting Ruelle pressure follows from (|2T| ) and (^) as 

\ 2 



as obtained in 



^ W (/3,P) 

ro|. 



pln(a + b) 




71 



L + a/ (bp) 



0(L" 4 ) 



(33) 



3.2 PBC — Configurations with a void 

An analysis of the eigenvalue A(/3) for a configuration with a single void 
of width R (i.e. R — 1 empty sites) in the PBC and ABC case (see figure 
|T]b) provides essential insights for dealing with more complex distribution of 
scatterers. In this subsection, PBC's are treated, in the next one ABC's. 

We start with a perturbative calculation for large R. Recall that the Haus- 
dorff dimension dn of the repeller is defined through ip(dn) = hiA((i#) = 0, 
and consider first (3 < dn, so that A(/3) > 1. The position of the left-most 
scatterer of the cluster is chosen as site number 1. The equations ([T^) for 
the scattering amplitudes read in this case 

A R U + (1) = aU+(N) + eb[U-(N)/e] 

AC/_(1) = 6*7+ (2) + aU-(2) 

AU + {2) = ea[U + (l)/e]+bU-(l) 

AE7_ (2) = W7+(3) + aU-(3) 

(34) 

ALL (A" — 1) = bU+(N) + ea [U-(N)/e] 

AU + (N) = aU + (N- 1) + bU-(N- 1) 

A R U^N) = eb[U+(l)/e] +aC/_(l) 

where the e(l/e) factor has been introduced to clearly display the structure 
of the following perturbative calculation. For a large void of length R the 
solution of (j34j) for a configuration in a closed system of length L = N + R — 1 
is expected to look like that for an open system of length N with N scatterers, 
studied in subsection |3.1| . This can be inferred from the first and last equation 
in fl34), where e = A 1_R is a small quantity for large R : the eigenvector 



components U + (l) and U^(N), corresponding to "entering" velocities, are 
expected to be linear in e. In a perturbation expansion in powers of e, 
the components U+(l) and U-(N) are vanishing to dominant order in e, 
corresponding to absorbing boundary conditions. The first and last equations 
in (^) can be considered as new boundary conditions replacing (|TTJ), and the 
remaining 2(N — 1) equations can be written in matrix form as 

AU = W°U + eA (U/e) (35) 
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where W°(k,i\k',j) is the matrix representation of ([13]) with = Rk-i = 1 
for the ABC case. Inspection of (S3) shows that the perturbation matrix A 
has the form 

A(k,i\V,j) = a5(k,2)5(i,+)6(k',l)6(j,+) 

+a5(k, N - l)5(i, -)6{k', N)5(j, -), (36) 



where S(k, I) is a Kronecker delta. The eigenvalue equation fl35|) can be 
solved by a perturbation expansion around the solutions {A , U%} of the 
open system, discussed in subsection ^Ij i.e. A = A G + eAi + • • • and Ui = 
U° + t\J\ + • • •, yielding the equations of 0(1) and 0(e) : 

(A -W°)U° = 
(A - W^U 1 + Axf/ = AU\ (37) 

whereas the new boundary conditions follow from the first and last line of 



A C/i(l) = aU° + (N) 

A U^(N) = aU°_{\). (38) 

Here A Q and U°(k) are explicitly given in fl32| ) and (]28f) with A = B = 
1. To solve the 0(e)-equations in Q3"7p, we need the left eigenvector V°(k) 
defined through V°W° = A V°. Using the symmetries in the explicit form 
of the 2(N — l)-dimensional matrix W°, it is straightforward to relate the 
components of the left and right eigenvectors, V±(k) and U±{k) respectively, 
with the result : 

V°{k) = U°(k- 1) = sin q (N - k + 1) 

V°(k) = Ul(k+1) = sinq k } (39) 

where q a = 7r/ (N + a/b) as given in (|31|). The eigenvalue in first order follows 
from (|37D by taking the inner product of (|3^) with yielding 



eAi = e i?^ K i (?) 2 (a + bf ' Rql (N ^ X) (40) 

where inner products are defined by 

(X\Y) = £ E A A j(/c)F i (/c). (41) 

fc i=±l 

Here numerator and denominator have been calculated from (^) , (|3T)| ) , (P5[ ) . 
and (|39"D with the result 

(^lAlf/ 1 ) = a^(2)E7i(l) + ayf(JV-l)E/i(JV) 

2a 2 [sing (iV - 1)] 2 /A ~ 2(a 2 /b 2 )(a + b)q 2 
(V°\U°) = cot q Q sin q Q N - N cos q a N 

~ N + a/b = n/q . (42) 

We conclude that the largest eigenvalue A(/3) and corresponding Ruelle pres- 
sure = lnA(/3) in a PBC — configuration with N scatterers and a void 
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Figure 2: Eigenvalue A(/3) at (3 — 0, for a configuration with only one void 
containing R — 1 empty sites and PBC, as a function of R. The upper 
(lower) curves correspond to the numerical results for N = 20 (N = 10); 
symbols (circles) give the prediction based on equations $4$)-(\m); horizontal 
lines indicate the exact eigenvalue A (j3 = 0) for a PRW on a lattice of size 
L' = N with ABC. For (3 = 0, £ = l/ln2 ~ 1.4 is the correlation length. 



of length R decays exponentially with a correlation length £ = 1/ ln(a + b) 
to the eigenvalue of a solid cluster of N scatterers with ABC's. 

In appendix A, we present a detailed calculation, exact for all R (R = 
1, 2, ■ ■ •). It yields for the largest eigenvalue 



A(q) = (a + &){!- {a/2b)q 2 + ■■■} 



where 



q = qi = ix / { N + 



1 + {a + bf- R 
\-{a + b) l - R 



(43) 



(44) 



A plot of A(gi) 
e = (a + b) l - R 



is shown in Fig. ^. The e-expansion of this result, with 
agrees with the perturbation result (fiTf). For R = 1 (no 
empty sites), the resulting wave number reduces to q\ = 0, as it should for a 
closed system. 

Again, equation (f44"D shows that q and A(g) decay within a correlation 
length £ = 1/ ln(a + b) towards the corresponding values q Q and A (q ) of an 
open system with a solid cluster of scatterers. 

allow us to calculate all chaos properties of the 



Equations P3| ) and 
configuration in figure (|l|,b) with the help of equation (21). 



3.3 ABC — configurations with a void 

We consider an open system with two solid clusters (see figure (|l],c)) with 
respectively M and M = N — M scatterers, separated by a void of length R 
and study the eigenvalue problem for (3 < d H so that A(/3) > 1. We expect 
that for sufficiently large R the M- and M-blocks become independent, and 
the components f/_(M) and U + (M + 1) of the eigenvector, corresponding 
to "entering" velocities, are small and can be treated as new ABC's in a 
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perturbation calculation. We therefore write the eigenvalue equation (p~3|) 
as a set of 2(M + M — 4) coupled equations for the scattering amplitudes 
{Ui(k)\77i{k)} with {k,i) E {l-,2+,2-,---,(M-l)-,M+|(M + l)-,(M + 
2)+, • • • (JV - 1) + , (N - I)", N + } of the general form : 





(45) 



with "absorbing" boundary conditions for the M- and M-blocks in the form 



UJ1] 



A R U„{M) 



A H U + (M + 1) 





eb [Z7 + (M + l)/e] + aU^M + 1) 

aU + (M) + eb [t/L(M)/e] 
0. 



(46) 



The block matrices W°(M) and W°(M) refer respectively to the M- and 
M-clusters, and have the same form as W°(N) for the iV-cluster in equation 
(PSj). The matrix A connects the block matrices to the "entering" states 
C7_(M) and Z7+(M+ 1) : 



A(ki\k'j) = a5{k,M -l)6(i,-)6(k',M)6(j,-) 

+a5(k, M + 2)5(i, +)5(k', M + l)5(j, +). 



(47) 



The boundary conditions ( fl"6|) couple the two blocks. These boundary equa- 
tions show that U-(M) and M + (M + 1) for large R can indeed be treated 
as small quantities of order e = A l ~ R , as in subsection |3]2. 

The analysis of this problem is similar to that in subsection |3.2| . The 
eigenvalue equation 





15]) can be solved by an expansion in powers of e. 
The eigenvalue problem to zeroth order in e reduces to two decoupled 
equations for the two isolated M— and M-clusters with ABC, reading 



(48) 



For sufficiently large M and M the largest eigenvalues of the block matrices 
are A (g ) and A (q ) (see equation fl43|) ). These are also eigenvalues for the 
whole system (|58|) with right and left eigenvectors : 

and {V°{q ) | 0} with q = ir/[M + a/b] 
and {0|F°(g o )} with q Q = n/[M + a/b}. (49) 



{U°(q_o) | 0} 
{0 | U°(q )} 



The right and left eigenvectors in (^) are again given by (pS| ) and (|3^) with 
iV replaced by M and M respectively. If M > M, then A(g ) > A(g ) 
and A G = A(q ) is the largest eigenvalue with the corresponding eigenvectors 
{U°{q ) | U°{q ) = 0} and {V°{q ) \ V°(q ) =0}. 

To linear order in e, the boundary conditions ( PS|) require for the compo- 
nents {f/ 1 ^) | F 1 ^)} : 



C/i(M) 
7J + (M+ 1) 

ZJi(iV) 





aU°_(M+ 1)/ A 

aU° + (M)/A 

0. 



(50) 
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Figure 3: Eigenvalue A(/3) at (3 = 0, for a configuration with two clusters of 
sizes 50 and 30, separated by a void of size R and ABC. The horizontal line 
indicates the exact eigenvalue A (j3 = 0) = 2 for a PRW on a lattice of size 
L = 50 with ABC; and = 0) = l/[21n2] ~ 0.7 is a theoretical estimate 
for the correlation length, in good agreement with the numerical results. 

The last equality on the second line follows as all components of U°(q ) are 
vanishing. First order perturbation theory for the largest eigenvalue yields 

eAi = e(V°|A|f/ 1 ) = 0, (51) 

as a consequence of (|50D . 

Second order perturbation theory yields a no n- vanishing result, propor- 
tional to e 2 ~ (a + b) 2 ~ 2R . Therefore the largest eigenvalue for the con- 
figuration of figure (|l],c) as a function of the width of the void R has the 
form 

A OR) ~ A + Const x e" {fl - 1)/? (52) 

with a correlation length £ = [21n(a + &)] _1 , a factor 2 smaller than in the 
previous subsection. 

In an ABC — configuration with two clusters containing M and M scat- 
terers respectively, and separated by a distance R, the largest eigenvalue A 
approaches the eigenvalue A G = A(q ) given by fl43"l) -(|44l) at an exponential 
rate. This A(q ) is solely determined by the largest cluster. The plot in Fig. 
[| shows the numerical solution of the eigenvalue problem. 

The main conclusions of the previous subsections, referring to (3 < dn, 
can be directly generalized to configurations with more clusters : 

(i) The largest eigenvalue A(/3) in a PBC — configuration (closed system) 
with at least one void larger than £ is equal to the largest eigenvalue for the 
corresponding ABC — configuration (open system). 

(ii) In any ABC — or PBC — configuration with clusters separated by dis- 
tances Rk > £ (k — 0, 1, • • • , N) the largest eigenvalue is A(q ) given by (f43"l ) 
with q Q = 7i/[M max + a/b], where M max is the number of scatterers in the 
largest solid cluster. 

In fact, we will use this last case (ii) to illustrate the localization process 
mentioned in the introduction. In section ^ we have seen that the largest 
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Figure 4: Invariant vector 7r + (r) = u + (r)v + (r) . The location of the two 
clusters has been indicated with dashed lines. 

eigenvalue A(/3) of the matrix w(/3) was dominating the dynamic partition 
function Z(j3,t\x Q ) at long times, and thus the Ruelle pressure (equations 
fll9|)-(|20D). On the other hand, we just found that, for (3 < da, A(/3) is 
determined by the largest cluster of the configuration, as if this cluster was 
isolated and surrounded by absorbing boundaries. It means that the trajec- 
tories which dominate in sum (|15|) are those which always remain localized 
inside the largest cluster. The other trajectories traveling throughout the 
system may as well be omitted. Thus the Ruelle pressure in large systems 
will not reflect the global structure of the system but only characterize the 
largest cluster present in the configuration. It is precisely this phenomenon 
that will be referred to as localization. This can be illustrated by calculating 
the invariant vector (see figure Only the states corresponding to posi- 
tive velocities are plotted here. States with negative velocities give the same 
result, up to a translation : 

7T-(r){u\v) = w_(r)t>_(r) = v + (r + l)u + {r + 1) = n + (r)(u\v) . (53) 

After an infinite time, all the probability of finding the particle is concen- 
trated on the largest eigenvector. 

To get some insight in how the Ruelle pressure varies with the con- 
figuration of scatterers, we may compare the largest eigenvalues A(/3) ob- 
tained for each of the three configurations of figure [|, keeping a fixed density 
p = N/L = (M + M)/L. The eigenvalues are respectively 



As p < 1 and (3 < du, i-e. a + b > 1, it is straightforward to show that 




(54) 



A a < A c < A fe . 



(55) 
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The largest eigenvalue is obtained when all scatterers are packed together in 
one solid cluster, while the smallest corresponds to the mean field configu- 
ration, where no cluster is formed and thus localization is not possible. In 
section it will appear that the largest eigenvalue of any other configuration 
falls between A a and We will show that, for large L, most configurations 
contain a largest cluster which will entirely determine the Ruelle pressure. 
Indeed, as we start to see with configuration (c), localization is not specific 
to very special configurations, but occurs more generally for most configura- 
tions. 

Finally, we stress that, in this section, only the case (3 < du was consid- 
ered. For (3 » 1, a complementary phenomenon occurs, i.e. localization in 
the largest void instead of largest cluster [|lj . 



4 Hausdorff dimension 



The eigenvalue equation in the representation (|T3| ) and the result of section 
|3.1| enables us to carry out an exact calculation of the Hausdorff dimension 
ds for an open LLG, which is defined as the root of ip{(3 = du) = 0, or 
equivalently as the root of A(/3 = dn) = 1 on account of (^T|). For a closed 
LLG there is no fractal repeller and du = 1. 

The important observation is that du is independent of the quenched 
disorder, and depends only on the total number of scatterers. This can be 
seen by combining equation (|13|) with the requirement A(/3) = 1. The random 
variables {Rk} disappear from the equation, so that du is the same as for 
the PRW in subsection |3.1| . It can be calculated by setting the right hand 
side of (521) equal to unity, and solving for (3. For large N (i.e. small q), the 



root is 




(3 = d H = 1- 



!"T I T I 7777^ ) +°(^ 4 ) (56) 

■^■closed 

where D = p/2pq is the exact diffusion coefficient of the one-dimensional 
LLG [II], and X c iosed — P^ — — p(plnp + glng) is the exact Lyapunov 



exponent for a closed LLG, as obtained in [10, Q 



5 Numerical method 

The remaining part of this paper describes the numerical diagnostics, in 
which numerical and analytical results will be compared. In this section 
we start with a description of the numerical method, used to calculate the 
largest eigenvalue of the large random matrix w(j3;x\y) in equation (H) for 
a fixed configuration of scatterers characterized by a certain system size L 
and number of scatterers N. Then the Ruelle pressure is obtained as the 
logarithm of this eigenvalue (equation (|2T|)). 

In one dimension, a recurrence formula can be found which allows us 
to compute numerically the exact value of the determinant of w{(3) — Al, 
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where 1 is the identity matrix. Then A can be determined as the largest 
root of the equation det\w{(3) — Al| = Or, using Newton's method. The 
recurrence formula is derived in appendix B. This method can be applied 
provided that L is not too large (less than 400). Indeed, if the system size is 
larger, numerical overflow problems occur. 

For large system sizes (L > 400), and for (3 1, the calculation of the 
determinant involves very large numbers that cannot be handled by worksta- 
tions. Under such circumstances, A has been determined by using Arnoldi's 
algorithm, which is an iterative method akin to Lanczos algorithm ||13|. Let 



w be an n x n matrix whose largest eigenvalue has to be determined. The 
idea is to scan rapidly the eigenvector space, and find a subspace U contain- 
ing the m most significant eigenvectors. Then we compute an m x m matrix 
if as a kind of projection of w onto the subspace U. The largest eigenvalue 
Ah of H associated with the eigenvector uh yields an approximation for the 
largest eigenvalue of w and Uuu for the corresponding eigenvector. If the 
result is not satisfactory, the whole process is repeated, taking Uuh as an 
initial guess. The method is explained with more details in Appendix C. The 
size m of the basis has to be tuned in order to optimize the efficiency of the 
method. 



6 Random configurations. 

In this section, we will illustrate that localization occurs not only in the 
special configurations considered in section 0, but in the majority of random 
configurations realized in large systems. 

To show this we generate some random configurations of a lattice of L = 
100 sites with a filling fraction p = 0.3, as shown in Figure |5j, and calculate 
the largest eigenvalue A(/3) and corresponding Ruelle pressure ip(/3) = In A(/3) 
numerically. We further determine the size M of the largest cluster in each 
configuration, and calculate the largest eigenvalue A(/3; g^) for an isolated 
M-cluster, using equations (f43|)-([44|) with q M = ir/[M + a/b}. The results 
are displayed in Table |I[ For comparison, the mean-field value for the system 
size L and number of scatterers N = pL is calculated from equation fl32|) and 
yields for (3 = 0, ^mf(O) = 0.2078 for ABC, and Vmf(0) = 0.2079 for PBC 
(where q = in equation (|2"7|). 

The Ruelle pressure lnA(/3) is in fairly good agreement with the estimate 
lnA(f3;q M ) for the three first configurations (a)-(c) (see Table [!]), i.e. as 
soon as the largest cluster size M is on the order of 5 sites. In other words, 
the dynamic partition function ( |i~5l) and Ruelle pressure ([2^) calculated from 
the subset of trajectories fl t that remain on the largest cluster for all time, 
give already a fair approximation to the actual Ruelle pressure defined by 
summing in ([15]) and ( |20| ) over all trajectories fl t that stay inside the domain 
T> for all time. This means that there is already a large degree of localization 
in configurations (a)-(c) and, to a lesser extent, in (d) — but not in (e) — , 
in spite of the small system size L = 100 and number N = pL = 30 of 
scatterers. 

On the other hand, if there is no large cluster (as in configuration (e), 
where M < 2), then A(/3;gM) is not a good estimate at all. The mean 
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Figure 5: Random configurations obtained for L = 100 and p = 0.3. The 
largest cluster is of size a) 6, b) 5, c) 5, d) 4, e) 2. The shaded areas represent 
the location of the scatterers on the lattice. 



configuration 


Ruelle pressure 


M 


estimate 


(a) 


0.59532 


6 


0.589 


(b) 


0.56275 


5 


0.549 


(c) 


0.55034 


5 


0.549 


(d) 


0.50916 


4 


0.481 


(e) 


0.35389 


2 


0.000 



Table 1: For the random configurations of figure^, with L = 100 and p = 0.3, 
we compare the actual Ruelle pressure ip(/3) with the estimate lnA(/3;gjvf) 
based on the largest cluster size M, for (3 = 0. PBC are used for configu- 
rations (a), (b), (c), and ABC for (d), (e). The corresponding mean-field 
values, ^mf(0) = 0.2079 (PBC) and ij> MF {o) = 0.2078 (ABC), do not pro- 
vide a sensible estimate. 

field value is slightly better but still is a poor estimate, as there are large 
fluctuations in the distances between the scatterers. 

6.1 Bounds on Ruelle pressure. 

In II 10, it has been shown that the Ruelle pressure is bounded by 
lnA(g M )~ln(a + &)-^7r 2 /M 2 < < l n (a + b) (fi < 1) 

/3Qnq)l/M < MP,P) < 

(57) 

where M is the size of the largest cluster and M of the largest void. The 
upper bound is a direct consequence from the inequalities 

a + b > 1 if < 1 (58) 
a + b < 1 if p > 1, (59) 

and we refer to || []J for more details. For (3 < 1 (respectively (3 > 1) 
the lower bound is the Ruelle pressure obtained by keeping only the largest 
cluster of the configuration (or respectively the largest void, bordered by 2 
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Figure 6: Among 1000 random configurations generated for L = 50, N = 25, 
we select the configuration corresponding to the largest (a) and smallest (b) 
A-value, as determined numerically. The shaded areas correspond to the 
location of scatterers. 



:+ 0.3 



> 0.2 




Figure 7: Invariant measure tt(x) of Eq. $25] ) for configuration (a) of figure 
|^. A thin line indicates the profile for the location of scatterers. 



scatterers) and using ABC. In subsections $.2[ |Q| , and || we have found that 
for (3 < 1, this is not only a lower bound but also a good estimate for the 
Ruelle pressure, as soon as L is large enough. A consequence is that, among 
all possible configurations, the configuration with all scatterers solidly packed 
in a single cluster gives the maximum value of the Ruelle pressure. On the 
other hand, the minimum value is obtained for the mean field configuration 
with scatterers spaced at regular intervals of length R = 1/p. 

As a confirmation, we have verified (see fig. that in a given set of 
1000 configurations, the largest eigenvalue for (3 < 1 did correspond to a 
configuration where all scatterers are essentially packed together, whereas 
the lowest value was obtained for an "almost mean-field" configuration, i.e. 
the distance between scatterers is more or less constant. This is in agreement 
with our expectations [|TJ. To illustrate localization in these configurations, 
we have plotted the invariant vector v + (r)u + (r) as a function of r (figures 
[7| and §). For configuration (a), the eigenvector is entirely localized on the 
large cluster on the left. Configuration (b), which is more of mean-field type, 
correspond more or less to an extended state. However, there is still partial 
localization in regions of higher than average density. 



18 



10 20 30 40 50 

r 

Figure 8: Invariant vector for configuration (b) of figure A thin line 
indicates the profile for the location of scatterers. 

7 Distribution of largest cluster size 

Up till now we have studied single configurations. In the next sections we will 
present results averaged over the disorder and compare them with the upper 
and lower bounds in Eq. fl5T|). To do so, we need to average equation (|57|) 
over all possible configurations of scatterers with a fixed density p (or a fixed 
N) and a fixed system size L. In order to determine (1/M 2 ) P and (1/M) p , we 
have used three different estimates for the distribution of the largest cluster 
size. Note that the distribution for the size of the largest void is the same, 
upon exchanging scattering sites and empty sites. The first estimate, which 
is the most crude one, is just the distribution for having at least one cluster 
of size M given by 

/ L-M -2 \ 

[ N -M j 
Pi(M) = L— , (60) 




This expression is valid only for large M and for periodic boundary condi- 
tions. The numerator in (^0]) represents the number of ways one can dis- 
tribute the N — M remaining scatterers among the L — M — 2 remaining 
empty sites, once a cluster of size M limited by two empty border sites has 
been put in one of the L possible locations. The denominator is the total 
number of possible configurations. 

The calculation of the distribution of largest cluster sizes can be improved 
in the following way: Let A(M) be the fraction of realizations with no cluster 
larger than M. Then 

A(M-l) = A{M) x [l-Pi(M)]. (61) 

Note that 1 — P\(M) is the probability that there is no cluster of size M. 
This recursion relation can be solved by iteration starting at M = N, where 
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Figure 9: Probability distribution for the largest cluster size M in a system 
of size L = 100, and p = 0.4. The dotted, dashed, and solid lines correspond 
respectively to P±(M), P 2 (M), and the direct measurement Ps(M) obtained 
by generating 20 000 random configurations. 

A(N) = 1. The result is : 

A(M) = n£ =M+1 (1 - P x (m)) . (62) 

The probability that the largest cluster size is M is then 

P 2 (M) = A(M) - A(M - 1) = Pi(M)II^ =M+1 (1 - P x {m)) . (63) 

This expression for P 2 (M) can be calculated numerically for each system size, 
by using the factorial expression given above for P\(M). This is the second 
expression which was used for the numerical evaluation of lower bounds. 

The third distribution P 3 (M) for the largest cluster considered here was 
obtained by simple generating a large number of configurations and finding 
the largest cluster for each of them. 

Figure ||| compares these three estimates for the distribution when p = 0.4, 
and L = 100. As Pi(M) is not bounded, it is cut off such that the probability 
is normalized. The distribution P 2 (M) has also to be cut off, otherwise it 
oscillates between unphysical positive and negative values at small M (but we 
stress again that the formulas ( pD| ) and (^) are not valid for small M-values). 
The third distribution P^(M) has been averaged over 20000 configurations. 

We now verify that the largest cluster size grows as InL for large L. 
Figure [1(^ shows P 2 (M) for increasing system sizes. The system size is varied 
from 100 to 10 8 and is multiplied by 10 between each successive estimation. 
We check that each time the size L is multiplied by 10, the maximum of the 
distribution is shifted to the right by a constant value. Another verification 
is made by plotting the second moment (1/M 2 ) 2 calculated with the P 2 - 
distribution as a function of (log 10 L)~ 2 (see figure |TT1) . 

As already discussed in 0, it is possible in the one-dimensional case for 
(3 < 1 and L not too large, that the Ruelle pressure is not determined by 
the largest cluster but by a "dominant" cluster with average density p + Ap 
intermediate between p and 1. This is indeed what is observed in numerical 
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Figure 10: Probability distribution P2(M) for the largest cluster size M for 
p = OA, (5 = 0, PBC, and for system sizes L increasing geometrically from 
100 to 10 8 . 
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Figure 11: Scaling properties of the second moment (1/M 2 ) 2 , calculated from 
figure IJZ], as a function of l/(log 10 L) 2 for L = 100 to 10 s , at p = 0.4 and 

p = o. 
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Figure 12: Probability distribution for the largest (solid line) and dominant 
(dashed line) cluster size in a system of size L = 100, for p = 0.4 and (3 = 0. 
These distributions were obtained by generating 20 000 random configura- 
tions. 



simulations. At (3 = and for L = 100, we have taken all segments of 
all lengths, measured the average density on each of them, calculated the 
corresponding Ruelle pressure using the mean field expression (|33D , and kept 
the one which gives the largest value. We call it the dominant cluster. For 
some configurations it coincides with the largest cluster, but not always. 
Figure [12] compares the distributions for dominant and largest clusters. They 
are different, the first one being slightly shifted towards larger values. Note 
that a cluster, dominant for a given (3, may not be dominant for another (3 
value. 

Another illustration of partial localization in high density regions instead 
of solid clusters was given by figure EL 



8 Average Ruelle pressure. 

Now we can use any of the three estimates for the distribution to calculate 
the average ( |57|) over all possible configurations and compare the resulting 
lower and upper bounds with the numerical measurements of the Ruelle pres- 
sure. For fixed L and p, a large number of scatterer configurations has been 
generated, and for each of them the largest eigenvalue has been calculated. 
The average of its logarithm yields a numerical value for the Ruelle pressure. 

We first consider the case where (3 < 1. We have collected data for 
(3 = 0, for which the Ruelle pressure equals the topological entropy. Figure 
P~3] confirms that numerical data are between the upper and lower bounds. 
The lower bounds based on any of these three distributions are qualitatively 
the same. 

The estimate for the largest cluster distribution, especially the one based 
on Pi(M), may seem rather crude. In fact, most of the system sizes that 
we are able to explore numerically are too small for the Ruelle pressure to 
be entirely dominated by the largest cluster, and no improvement of the 
estimates for P(M) is likely to make the quantitative agreement better. 
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Figure 13: Ruelle pressure if) (ft = 0) (solid line) for L = 100, as a function 
of the density p, compared with upper (dotted line) and lower bounds (dashed 
lines). For comparison, the mean field prediction has also been indicated 
(dashed-dotted line). 

However, when the system size is increased (L = 10000), we check -within 
the precision of our measurements- that the lower bound based on the largest 
cluster approximation is indeed a good estimate of the Ruelle pressure (figure 



Hp . It should be noted that at ft = 0, the Ruelle pressure is independent of 
p and q, and thus these results are valid both for large or small q. 

One can also verify that for a given system size, the mean field prediction 
gives a value much lower than the average, whereas the configuration with 
all scatterers are solidly packed together has a Ruelle pressure almost equal 
to the upper bound. This is illustrated in figure [15| for L = 400. Thus for a 
configuration with L = 400 and a given density, any value between the two 
dashed-dotted lines of figure [T^ can be realized. 

Now we address the case ft > 1. Tables || and ^| give some numerical 
values for the measured or estimated Ruelle pressure, and for the mean field 
prediction (|33|), in the case of L = 100 and ft = 2. When q is large (table 0), 
the lower bound is not only a lower bound but also a good estimate of the 
Ruelle pressure, and is much lower than the mean-field value. The numerical 
data are also displayed in figure |16[ 

When q is small (table [| figure [H]), the situation is reversed. A good 
estimate is obtained by using the mean field theory result, while the theo- 
retical lower bound differs significantly from the measured Ruelle pressure. 
We may roughly estimate under which condition the lower bound will be a 
better estimate than the mean field theory by comparing the corresponding 



eigenvalues in ( p7|) and (p2|) with R — l/p : 

q P/M < ( Q + b y ( 64 ) 

or with b = q 13 

M > -j^r-, (65) 
p\ ln(a + b) | 

If p > q, then for ft » 1 we have a » b. Localization will dominate over 
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Figure 14: Ruelle pressure ip((3 = 0) (solid line) for L = 10000, as a function 
of the density p, compared with upper (dotted line) and lower bounds (dashed 
lines). For comparison, the mean field prediction has also been indicated 
(dashed-dotted line). 




Figure 15: Ruelle pressure ip((3 = 0) as a function of the density p for L = 
400. The solid line is an average over 10000 configurations. The dashed- 
dotted lines are the extreme values obtained for specific configurations, namely 
the 'mean field' configuration and the one where all scatterers are packed 
together in a solid cluster. The upper bound ip(0) = In 2 is also indicated. 
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Figure 16: Ruelle pressure as a function of the system size in the case of 
strong backs cattering (p = 0.2 ) for p = 0.2 (o) or p = 0.8 (n ). Lower bounds 
(solid lines) were obtained using a numerical determination of the largest 
cluster size distribution. 



P 

Numerical value 
Lower bound 
Mean field 



0.2 0.5 0.8 

0.02809 0.06583 0.12211 

0.02851 0.06682 0.12406 

0.07713 0.19283 0.30853 



Table 2: For p = 0.2 (large backs cattering probability), the Ruelle pressure, 
evaluated numerically by averaging over 10 000 configurations, is compared to 
the lower bound and the mean field value, for different densities of scatterers 
and a system size L = 100. 



mean field if 



M > 



ln&j 



In 9 1 



(66) 



p|lna| p| lnp| ' 

i.e. the system size L must be typically larger than exp(l/p lng/lnp). On 
the other hand, at fixed L, (3 must be larger than 



P>1 + 



(lng)|ln(l -p)| 
(p In p + q In q) pM 



(67) 



where M ~ InL. This results from an expansion in powers of /3 — 1 [|1]. It 
appears that, if q is small, it is necessary to go to very large L and/or (3 
values to see the occurrence of localization. In the case of figure [T7|, (3 = 2.0 
and L = 400 are not large enough to see this phenomenon. Figure [18] shows 
that this is also the case for a different density p. 

The picture of what happens when p > q can still be refined. The 
crossover between trajectories extended over the whole lattice (mean field) 
and trajectories localized in the largest void (lower bound) is in fact not 
direct. Some intermediate semi-localization may occur. As mentioned in 
section |^, this was already evidenced for one-dimensional systems at (3 < 1 
(and for all p, q values) |L). We will show now that a similar phenomenon 
occurs for (3 > 1 when p > q. We then have not only a » b, but also 
a « 1 (j3 > 1). Thus (i) on the one hand the particle has a tendency to 
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Figure 17: Ruelle pressure as a function of the system size in the case of 
weak backs cattering (p = 0.8) at a high density p = 0.8 (A). The lower 
bound (dashed line) was obtained using a numerical determination of the 
largest cluster size distribution. The mean field prediction (dashed-dotted 
line) is a much better estimate for the Ruelle pressure. 
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Figure 18: Ruelle pressure as a function of the system size in the case of weak 
backs cattering probability (p = 0.8) at a low density p = 0.2 (A). Again, the 
mean field prediction (dashed-dotted line) is a much better estimate for the 
Ruelle pressure. 



p 0.2 0.5 0.8 

Numerical value 0.07218 0.17361 0.29033 

Lower bound 0.20584 0.48190 0.91901 

Mean field 0.07713 0.19283 0.30853 

Table 3: For p = 0.8 (weak backs cattering probability), the Ruelle pressure, 
evaluated numerically by averaging over 10 000 configurations, is compared to 
the lower bound and the mean field value, for different densities of scatterers 
and a system size L = 100. 
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Figure 19: Special configuration, in which a trajectory exploring the whole 
low density region may have a higher weight if p > q than a trajectory 
remaining in one of the voids. 



escape from the void because back-scattering is weak; (ii) on the other hand, 
as a << 1, the free propagation is till favored over forward propagation 
through a scatterer. If a ~ 1, the second effect is negligeable and the mean 



field prediction will be appropriate, as it is the case in figures |17| and [L8| If 
a « 1, the competition between the two effects may eventually promote 
some intermediate configurations. As an illustration, we will consider the 
configuration of figure [19] : the whole lattice is solidly filled with scatterers 
except in the low density region of size R. In this region, n isolated scatterers 
are placed at equal distance from each other. We will compare the weights 
of a trajectory 7\ undergoing only free propagation and back-scattering, i.e. 
confined in a void in a strict sense between two isolated scatterers; and a 
trajectory T 2 going through the n isolated scatterers and thus exploring the 
entire low density region of size R. 

During t time steps, T x will undergo t/[R/ (n + 1)] back-scatterings, while 
T 2 will undergo t/R back-scatterings and tn/R forward-scatterings. Thus 
T2 has a higher weight than T\ if 

at atn g t(n+l) 

qf-RpP r > qP^r- (68) 

or equivalently p > q. This shows that as soon as p > q, localization may 
occur not in the largest void but in a large low density region. However, 
when L increases, we speculate that it is more and more likely to find a large 
void that will nevertheless dominate the result. 



9 Conclusion 

We conclude this paper with a number of remarks: 

1) We have shown that the Ruelle pressure of Lorentz lattice gases in 
the limit of infinite systems is completely determined by rare fluctuations in 
the configuration of scatterers. Thus it carries no information on the global 
structure of the disorder. Numerical studies allowed us to show that this also 
holds for finite but large systems, for all but a small range of (3 values. Then 
we were able to predict quantitatively the Ruelle pressure for all f3 except for 
a small region around (3 = 1. 

As already discussed in ||, these results can be immediately generalized 
to a whole class of diffusive systems with static disorder. The special case of 
a continuous Lorentz gas has been briefly considered in 0. 

Localization phenomena in fact appear very often in physics, as soon 
as there is some competition between energetically favorable configurations 
and entropic effects. Some localization effects were already shown in the 
framework of the thermodynamic formalism for deterministic maps or 
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multifractals Qj. However, it is the first time that such effects are evidenced 
in hard spheres systems as resulting from the quenched disorder. Localization 
occurs on the most extreme fluctuation of the disorder. For infinitely large 
systems, this fluctuation may be arbitrarily large, which allows for a very 
pronounced effect. 

2) The Ruelle pressure is a characteristic of the dynamics of the system 
not only for the isolated (3 values where a direct interpretation can be given 
(see section §), but also as a whole. In the same way, in a power spectrum, 
only some of the points may receive an individual interpretation, but the 
whole structure of the spectrum is interesting. An open question is to know 
if enough information has been kept in the region around (3=1 where de- 
localization does not occur, in order to be able for example to reconstruct 
the structure of the disorder from it. To this respect, it would be interesting 
to rescale the region around (3 = 1 as the system size increases, in order to 
prevent it from shrinking to zero in the thermodynamic limit. The scaling 
of this region with the system size has been estimated in [|TJ. 

When localization occurs, the information contained in the Ruelle pres- 
sure concerns the properties of individual scatterers. More precisely, the 
knowledge of the (3 regions in which one given type of scatterers will dom- 
inate yields a measure of what could be called the isotropy of the different 
scatterers. 

3) Numerical studies have been performed only in one dimension. Ana- 
lytical results showed that the localization occurring in the thermodynamic 
limit and the extension of the delocalized region around (3 = 1 can be gen- 
eralized to higher dimensions. Localization occurs also for finite but large 
systems. The only difference with the one-dimensional case is that now the 
lower bound may not be a good estimate for the Ruelle pressure at finite size. 
Indeed, the lower bound chosen here was based on localization in hypercubic 
domains. In fact, it may happen on domains with much more general shapes. 
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Appendix A 



In the case of N scatterers packed in a single cluster in a system of size 
L = N+R — l with PBC, a perturbative calculation for the largest eigenvalue 
A has been presented in section |3]^. Here we calculate the largest eigenvalue 
directly from the set of equations (|3~4l). 

First notice that the bulk equations in (|3^) impose the relation (pE3|) be- 
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tween A and the wave number q, as can be found in a similar way as equation 
fl2"7D- Now the wavenumber q has to be determined from the boundary equa- 
tions 

A R U + (1) = aU+(N) + bU-(N) 
AU-(1) = bU + (2) + aU-(2) 
AU+(N) = aU + (N - 1) + bU-(N- 1) 
A R U_(N) = W7 + (l) +aC/_(l) (69) 

We search for a solution of the form 

L/+(Jfe) = Aexp[«g(A; - 1)] + c.c. 

U-(k)=B exp [«<jf(iV - fc)] + c.c. (70) 

By inserting this form into fl69"|) , we obtain four equations which determine A 
and B (complex numbers). Non trivial solutions exist only if the determinant 
of the system vanishes. This yields a transcendental equation for q which can 
be solved using the ansatz 

^n^TTS- < 71 > 

This ansatz is based on the assumption that the present case is similar to 
ABC, as soon as R is large enough. This will be correct only for (3 < 1. 

The determinant is expanded in powers of q. Replacing A by its expression 
(|43|) in terms of q, we find: 

Det = -4q 2 (5b-a + b)(A R + a + b) 

x ~8b(A R -a-b) - {a + b){A R + a - b)] . (72) 

Setting this determinant equal to zero yields two solutions. We select the 
one that gives the smallest q, i.e. the largest eigenvalue, and find expression 

m 



Appendix B : Recurrence formula for the determinant of w 



For a one-dimensional system, the transition matrix has a form such that 
an exact recursion relation can be found to calculate its determinant. We 
will first illustrate it for ABC. The transition matrix is of the form 



/ A h 
bi A 

a 2 



w 



«2 

A 

b 2 



h 
A 

«3 



a?, 

A 

h 



bi 
A 



A 

bi-2 



b L -i 
A 

a-L-i 



cll-i 

A 

b L -i 



\ 



b L 
A J 
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If the fc-th scatterer is located in site j, we define as the determinant 
of the 2j — 1 first lines and columns. An auxiliary quantity Ek is obtained 
from Dk by removing the last column and the last but one line. Then the 
following recursion relations hold : 



Dk+i 
Ek+i 



Using the initial conditions 



Ex 



A 2Rk D k - bE k+1 
bD k + a 2 E k . 



0. 



(74) 



(75) 



we can iterate the equations ([74]) and find the determinant for the whole 
matrix w as 

Det(w) = A 1+2Rn D n . (76) 

The recursive calculation is performed numerically. For (3 = (i.e. A ~ 
2), it is possible to calculate the determinant for matrices up to size 4(L — l) 2 , 
beyond which the method is spoiled by numerical overflow problems. 

In the case of PBC, the recursion relations are slightly more complicated. 
The matrix w is now of size 4L 2 and reads 



w 



( A 

A b 2 a 2 
a\ b\ A 

A b 3 a 3 
a 2 b 2 A 

A 6 4 
a 3 b 3 A 



a-L b L \ 



\ h ai 



A b L -i a L -i 
b L -2 A 

A b L a L 
a-L-i A 

A / 



(77) 



We choose the origin of the lattice such that there is no scatterer in sites 
1 and L. For PBC, this is always possible if there are 2 holes in a row 
somewhere in the configuration, which will be the case for all configurations 
if p < 0.5. Then b x = b L = 0. 

It can be shown that the determinant of w is given by 



Det(w) = A 2Rn+1 D n - bD* N - a 2 E* N + 2A(-A) 



L " 1 a 7V , 



(78) 



where Dn and En, D* n and E* N , are obtained from the recurrence formula 
( |T41 ) with respectively the initial conditions ([TJJ) and 



D * = A 2(Ro-i) b2 _ bE * 
El = b 2 2 -a 2 2 . 



(79) 
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Appendix C : Arnoldi's method 



Consider annxn matrix A. 

We start from a guess u Q for the right eigenvector. In a classical power 
method, one would apply the matrix A repeatedly to this vector until it is 
aligned with the eigenvector associated with the largest eigenvalue. With 
Arnoldi's method, we built a basis of m vectors that will span the vector- 
space more rapidly. If j vectors {u , u±, ■ • • , Uj-i} of the basis have been 
obtained already, Uj is defined as follows. First matrix A is applied : 

u'j = Auj^. (80) 

Then u'- is orthogonalized with respect to the j first vectors : 

3-1 



;. - £ [u' r u k ] u k . (81) 



k=0 



Finally, Uj is equal to the normalized u". This process is iterated until m 
vectors are obtained. 

A reduced m x m matrix H is defined by 

Kj = « i_1 u J ', (82) 

such that 

A = UHU T + B, (83) 

where U = {u Q , ui, • • • , u m }, and B is expected to be small. Notice that the 
definition of H implies that hj+ij is the norm of u", and hij — if i > j + 1. 
As a consequence it is straightforward to calculate the determinant of H, 
and thus to find its largest eigenvalue A by Newton's method (we know 
that it is smaller than a + b) and the associated eigenvector The first 
approximation for the largest eigenvalue of w is taken to be A associated 
with Uhu- If it is not satisfactory, the whole process is repeated, taking 
UtiH&s a new initial guess. The size m of the basis has to be tuned in order 
to optimize the efficiency of the method. 

A difficulty that this method shares with other iterative methods occurs 
when the largest eigenvalue is almost degenerate with the next smaller one. 
Then we may by mistake converge towards the second one. However, this is 
not a real problem as long as we are interested in the eigenvalue itself. 

It should also be noted that, as usual, convergence theorems exist only 
for symmetric matrices, whereas the method has been applied here to non- 
symmetric matrices. 
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